Nonlocal optical generation of spin and charge currents on the surface of magnetic 
insulators using total absorption and surface plasmons 
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We study the nonlocal spin and charge current generation in a finite metallic element on the 
surface of magnetic insulators such as yttrium iron garnet due to the absorption of the magnetic 
surface plasmon (MSP). Whereas a surface plasmon is completely reflected by a metal, an MSP can 
be absorbed due to the absence of backward states. The injection of MSP generates a voltage in 
the longitudinal direction parallel to the wave vector, with the voltage proportional to input power. 
If the metal is a ferromagnet, a spin current can also be induced in the longitudinal direction. Our 
results provide a way to improve upon integrated circuits of spintronics and spin wave logic devices. 
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Spintronics, which provides a promising way to over- 
come the limitation of charged-based electronics in data 
storage and processing, has received enormous attention 
recently A crucial process in spintronics is the injec- 
tion of spin and charge currents between structures, at 
least one of which is magnetic. The injection current is 
reduced when there is an impedance mismatch across the 
junction. Improving the injection efficiency is of the ut- 
most importance. The injection can be local as in vertical 
structures or nonlocal as in planar structures. Nonlocal 
structures hold the promise for large scale integrated cir- 
cuits applications. Current nonlocal spin transfer struc- 
tures [2| exploit diffusive spin currents. In this structure, 
the distance travelled by the signal is equal to the spin 
diffusion length l s f, and is of the order of 1000 A. Re- 
cently the focus of spin injection in vertical structures 
has shifted to structures involvingmagnetic insulators 
such as yttrium iron garnet (YIG) |3J , where the loss may 
be smaller. 

The exploitation of magnetic surface plasmons (MSPs) 
in spintronics have not been discussed previously but the 
MSPs offer many intriguing potentials. MSPs on bulk 
insulating magnetic surfaces are also known as Damon- 
Esbach [4| magnetostatic modes. They are coupled 
modes of the surface spin wave and the electromagnetic 
(EM) wave and are the magnetic analog of the better 
known (electric) surface plasmons. These states differ 
from bulk spin waves in that the frequency of these states 
are in the band gap region and corresponds to nonprop- 
agating bulk spin waves. They differ from surface plas- 
mons in that they are nonreciprocal, because of the bro- 
ken time reversal symmetry of magnetic systems. This 
leads to the absence of backscattering, which is also ex- 
hibited by electronic "edge" states in the quantized Hall 
effect [H and in topological insulators 0. Stimulated 
by recent focus on plasmonics, there has been interest in 
MSPs which can be generated on the surface of magnetic 
insulators such as YIG Q ■ 

When an EM wave falls on a metal, it is strongly re- 
flected due to the high impedance mismatch; very little 



signature of this wave is left in the metal. A surface plas- 
mon has an EM wave component and thus is completely 
reflected as well when it strikes a metal. In plasmonics, 
the injection of surface plasmons into metals has never 
been considered. Whereas a surface plasmon is nearly 
totally reflected, a MSP incident on a metal is absorbed 
now that the reflected surface wave is inhibited by the 
broken time reversal symmetry. 

In this Letter, stimulated by a desire to improve upon 
the spin injection efficiency, to reduce the loss of the sig- 
nal over distance of nonlocal structures, by recent inter- 
ests in spin wave logic devices Q and by experimental 
results of voltage generation in a metal film M on YIG 
under microwave excitation we study the injection 
of an MSP into M (which may or may not be ferromag- 
netic), as is illustrated in Fig. 1. We found that the 
MSP can be efficiently injected into M by the incom- 
ing electric and magnetic field regardless of the interface 
"impedance mismatch" , generating a static voltage due 
to the Hall effect in the metal and, also, a spin current if 
M is a ferromagnet. This provides for the electrical de- 
tection of the MSP. In addition, since the transport of the 
MSP is ballistic, the loss of the signal over propagating 
distance is much reduced. The distance travelled by the 
MSP is macroscopic, of the order of a centimeter. It is 
larger than the spin diffusion length of current nonlocal 
spin valve structures by five orders of magnitude. Finally, 
due to the partial spin wave nature of the excitation, the 
MSP can be easily coupled to spin wave logic devices and 
provides an efficient connection between different circuit 
elements. We now describe our results in detail. 

The nonlocal geometry we have in mind is illustrated 
in Fig. 1. Magnetic surface plasmons can be generated 
at the surface of a magnetic insulator with an external 
oscillating RF field, such as from a strip line, away from 
a finite metallic element M deposited on the surface of 
a magnetic insulator. In the following we shall illustrate 
our ideas with parameters such that the magnetic insu- 
lator is YIG. 

When fully magnetized, a magnetic material is de- 
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FIG. 1: (color online) Schematic diagram of the geometry for 
the injection of an MSP into a metallic element M. The MSP 
has wave vector k along x and the electric field E parallel to 
the magnetization M s . 
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FIG. 2: (color online) The dispersion of the MSP at the sur- 
face of YIG. The yellow (cyan) region denotes the frequency 
range where only MSP states with positive (negative) k x is 
allowed, while the MSP states with negative (positive) k x are 
inhibited, due to the the broken time reversal symmetry. 
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where fi and \j! are of the resonance form: /i = l+/ m (/o+ 
iaf)/[(h+iaf) 2 -f\ji' = -f m f/[(fo+iaf) 2 -f] with 
a charcteristic bare spin wave frequency /o = ^H . Here 
Ho is a sum of the external field H ext and the anisotropy 
field H a , 7 is the gyromagnetic ratio. f m = 7 4ttM s mea- 
sures the coupling strength of the magnetic material with 
the EM waves and a denotes the magnetic damping. For 
our calculation, we use Hq = 900 Oe, 4irM s = 1750 Oe, 
a = 7 x 10~ 3 , and the permittivity e s = 15 + i7 x 10~ 3 . 

We assume a geometry where the YIG surface and the 
film is in the x — z plane with the applied magnetic field 
and the magnetization M s along the z direction. We 
assume M to have straight edges parallel to the x and 
the z axes. The frequency of the MSP at the air- YIG 
interface as a function of k x for k z = for the above 
parameters is illustrated in Fig. 2. 

There are four frequencies that are of interest. The 
bulk spin wave gap at long wavelengths occurs between 
the spin wave frequency f sw = y/ /o(/o + fm) and the 
frequency f b = / + f m . The MSP frequency is in 
this band gap region where the bulk spin wave states 
are no n-propagating. The other tw o frequencies are 
fu = Vi e sfo{fo + fm) - / 2 ]/(e s - !)> where the disper- 
sion with negative k x starts and the frequency f s = 
fa + /m/2, at which the effective magnetic permeabil- 
ity ;U e ff = /i — jj! = — 1 and the dispersion with positive 
k x ends. For the above choice of parameters, they are 
given by f sw = 4.324 GHz, /„ = 4.425 GHz, f s = 4.97 
GHz and f b = 7.42 GHz. The MSP frequency starts at 
f sw . In the frequency ranges f sw < f < /„ (yellow re- 
gion) and f s <f<fb (blue region) the frequency of the 
MSP is defined only for one sign of the wave vector k x 
[loj . Because of the broken time reversal symmetry of 



magnetic systems at these frequency ranges there are no 
MSP states for wave vectors of the other sign. As can 
be seen from this figure, the group velocities Vmsp of the 
two branches of the MSP are of the order of 10 8 cm/s and 
the speed of light, respectively. The distance travelled is 
controlled by the damping coefficient a and is of the order 
vmsp / (ccw). For the two branches, they range from 10cm 
to 1000 cm. The y extent of the MSP, ((uj) 2 /c 2 ~kl)-^ 2 , 
is of the order of a millimeter. 

We next discuss the interaction of the MSP with the 
metallic film. The usual formulation of spin and charge 
injection focuses on the matching of the spin and charge 
current at an interface. In the present case there are 
no charge current; the boundary occurs in air and there 
are no discontinuity in the magnetic properties in the 
YIG substrate. A MSP is a mixture of a spin wave and 
an EM wave localised near the interface. The MSP is 
completely described by its EM field, the information 
of the spin wave is implicitly included in the magnetic 
susceptibility. The interaction of the MSP and M can 
thus be formulated in terms of its EM fields. 

When an EM field strikes a metal, usually it is strongly 
reflected. This is what happens with the electric surface 
plasmon. For the MSP, there are no states with negative 
k x s that it can be reflected to. We have considered the 
possibility that the MSP is stopped at the inteface, then 
goes along the edges of M and around it. This will involve 
states with finite k z and with k x — 0. We find that 
states with k z ^ are not possible. This can be seen 
from the following physical argument. For a nonzero k z , 
the transverse electric (TE) and the transverse magnetic 
(TM) modes are coupled; the EM wave in YIG have a 
finite component of the electric field E y in the direction 
perpendicular to the film. A finite E y induces electric 
charges at the interface and implies that the state will 
also have an electric surface plasmon character. This 
is not possible unless one of the dielectric constants is 
negative. 
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We have also considered the possibility of the wave 
going up at the air-M boundary. To this end we have 
studied states at the air-YIG interface localized near M 
with imaginary in plane (along x) wave vectors but real 
wave vectors in the direction perpendicular to the inter- 
face (along y). We found that there are no states going 
up from the air-YIG interface. The only states possible 
are those coming in from above. Thus there is no surface 
states for the incoming MSP to escape to and this state 
can be injected into M. In the following to illustrate our 
result we consider the case with the MSP wave vector per- 
pendicular to the direction of the magnetization (k z = 0). 
The electric field is then along z. The tangential electric 
and magnetic fields E z , H y and the perpendicular field 
B x inside M is the same as that of the incoming MSP 
as there are no reflected wave. We denote quantities in 
M and YIG with superscripts M and Y; the quantities 
in air have no supersctripts. The tangential electric field 
inside M, E z , generates a current with current density 
given by J^ 1 = uE^ 1 where a is the conductivity. Be- 
cause of the Hall effect, the magnetic field along y and 
this current along z generate an electric field along x of 
magnitude 

= -B%B?J?, (2) 

where Rfj w — 1/ (ne) (n is the density of the electrons) is 
the Hall coefficient of M, and the superscript H denotes 
the E field due to Hall effect. As we emphasize at the 
beginning of this paper, the fields only depend on the 
intrinsic properties of the metal and are not functions 
of the interface properties and any possible " impedance 
mismatch" . The product of E^ 1 and B^ and hence E x 
possess finite time independent components. In general 
M has no electrical leads attached to it and is in an open 
circuit situation. Charges will be induced at the surface 
until an electric field of magnitude (E x ) is developed and 
the current stops flowing. This corresponds to a voltage 
U x = dx(E^) where L x is the length of M in the x 
direction. The total voltage is controlled by the spatial 
dependence of the EM wave inside M. We show next that 
the electric field decays in M with a length scale of the 
order of an effective skin depth. 

The electric field in a medium is described by 
Maxwell's equation flE(r) = where Q = V x ( u _1 /eV x 
— u> 2 ; /i is the magnetic permeability ; e, the dielectric 
constant. Here we have used units so that the speed of 
light c=l. In principle, the electric field inside M/YIG 
can be calculated from its value at the boundary with a 
Green's function technique [Tlj from the equation 

E(r) = -JdS- [(//"VeV x E(r s )) x G(r, r s )* (3) 

-(n-'/eV xG(r,r s )y x E(r s )}. 

where the surface integral is carried out at the YIG-air 
YIG-M boundary. G is the Green's function defined by: 



n'G(r,r') = 5{r' - r). The operator f2' = Cl M for y > 0, 
il' = SlyiG for y < 0. G can be constructed from a linear 
combination of products of photonic states including the 
M-YIG MSP states as G = £ fex \n >< n|/(Ag - ul(k x )) 
where u) n , \n > are the eigenvalues and the eigenstates of 
il' . The Green's function is dominated by contributions 
from the MSP states at the air-M- YIG interface with 
energy Lo n (k x ) = k 2 when the denominator is zero. Thus 
the spatial dependence will be determined by these states 
which we study next. To illustrate the physics we first 
consider the simpler case where the thickness of M is 
larger than the spatial extent of the MSP. The MSP state 
at the air-YIG interface is coupled to the MSP state at 
the M-YIG interface which can be constructed as follows. 

MSP for M/YIG: The properties of the interface MSP 
state is obtained by matching ify and B± of the plane 
wave states of YIG and that in M. We obtain the disper- 
sion for the YIG-M MSP state given by 

k 2 xM a -uj 2 e M /[(b' - b + a? I a! - a']. (4) 

Here a = tf x (YIGf), b = ^ y (YIG); a' = n~} x {M) ,b' = 
/i~y(M). e, 6m, are the dielectric constant of YIG and 
M. These states can be excited by the incoming MSP and 
decays exponentially with a length scale of the order of 
inverse skin depth 1/k = l/Im[(wejv//a') 1 ^ 2 ] a way from 
the air/ YIG-M/ YIG boundary. 

MSP for a trilayer structure: So far we have assumed 
that the thickness of M is larger than the spatial extent 
d y of the MSP. d y is of the order of a mm. Current 
experiments are carried out with the thickness t of the 
metallic film of the order of nm and much less than d y . It 
is necessary to take the finite thickness of M into account. 

We have calculated the interface states in the 
air/M/YIG trilayer structure by matching the boundary 
EM fields at both the air/M and the M/YIG interfaces. 
We found that in the limit when t << 1/k, k x as given 
by Eq. (0| is modified so that e m is replaced by an effec- 
tive dielectric constant that is reduced by an amount of 
the order of tn. Furthermore, the real and the imaginary 
parts of k x are opposite in sign. In place of Eq. (jH) the 
dispersion is now given by 

k xM /[io 2 e M t} a -1/(1 + a + b); 1/(1 + a - b). 

The two solutions have opposite signs for the real part of 
k x . The voltage generated is determined by the imaginary 
part of the wave vector and is given by 

U x = Rh<?B y E z min{l, L x ) (5) 

where Z -1 ss 2oj 2 Im[eM]t/(l + a — b) is of the order tn 2 . 
For t « lOnm and a skin depth of the order of a 
micron, we obtain I « 0.1mm. We next estimate the 
magnitude of this effect. 

We take M to be Au with a Hall coefficient Rh = 7.25 x 
10~ n ?™ 3 /C, and a conductivity a = 4.55 x 10 7 /(ohm — 
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to). Thus E% = 3.3 x lQ- s B y E z . We estimate the prod- 
uct BE from the power P of the input coplanar waveg- 
uide. We assume that most of the power is expanded in 
generating the MSP. The energy current 



S = EB/hq = P/A 



(6) 



where A is the area of the wave. For current experi- 
mental systems, the width in the z direction is of the 
order of a mm; the y extent of the MSP is also of 
the order of a mm. Thus A rj 1mm 2 . We obtain 
< E z B y >« AttPc/ (IOvmsp)- Using the above estimate 
of I and for P of the order of a watt, U x is of the order 
of a micro- volt. 

Most metals are paramagnets. Because of YIG, or in 
the presence of an external magnetic field, the electrons 
in M has a small magnetization. A nonfcrromagnctic 
metal usually do not possess local moments and we do 
not expect exchange coupling between the YIG spins and 
this small magnetization in the paramagnet. The electric 
field E z generates an AC current inside M. As the elec- 
tric field oscillates from parallel to opposite the direction 
of the induced magnetization, a residual DC current can 
also generated if the resistances along the two directions 
are different. This can happen if there is a significant 
amount of spin orbit impurity scattering. The resistance 
along the two directions will then be different, resulting 
in a net average current / voltage in this transverse direc- 
tion. We next consider the situation when the metal is a 
ferromagnet. 

The Hall coefficient Rh of ferromagnets contains a 
component due to the anomalous Hall effect and the 
magnitude of the above effect can be modified. In ferro- 
magnets a spin current is also induced. At the air/YIG- 
M/YIG boundary, the perpendicular component of the 
magnetic field B x and the tangential component H y are 
continuous. Whereas the magnetization density in air 
is zero, there is a finite change in the x, y components 
of the magnetization density at the boundary inside the 
ferromagnet. We can obtain the magnetization density 
from the equation M(M) = [B(M)/fj, - H(M)]. M x<y 
act as sources to create magnetization currents in the fer- 
romagnet. From these we get the magnetization current 
given by Jm = i^M with 

M x = (a' - l/(io - b l2 /a')B x + ib'B y /{[i a'), (7) 

M v = {B y [l - 1/Guoa')] - ib'B x /a'} / Mo . (8) 

Associated with the change in M X}V there is a static 
change in the longitudinal magnetization M z to sec- 
ond order: AM, = 



Ml 



Ml - Ml 



M . This in- 
duces an additional charge current given by[l2| J' = 
-DmV(AM z ) where Dm is the diffusion constant that 
relates the charge current to the longitudinal magnetiza- 
tion change. In the open circuit situation, this provides 
for an additional contribution to the DC potential given 
by 
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FIG. 3: (color online) The output voltage from the spin 
current normalized by the incoming power as a function of 
frequency for the higher frequency branch of the MSP for 
external magnetic fields H ex t = 450Oe (dotted line) and 
H ext = 850Oe (solid line). 



We next estimate the magnitude of U as a function of 
frequency. 

The ratio Dm/& is of the order of l/{N(Ep)iiBe) 
where N(Ep) is the density of states at the Fermi sur- 
face; /zb, the Bohr magneton; e, the electric charge. 
From eq. © we get < B 2 >= P/j, /(vmSpA). We take 
M to be permalloy with a density of states N(Ep) « 
3 x 10 28 /ev/m 3 . For the spin wave parameters, we found 
recently fl3| that the fmr frequency of thin permalloy 
film as a function of frequency and external field H ext 
can be fitted by the typical spin wave formula f sw with 
satutaion magnetization (J-oMq = IT, and total field 
i?o = H ext — 126. 60e. Our estimates for the voltage U 
normalized by the incoming power P as a function of fre- 
quency for the upper branch of the MSP are shown in 
figEl for external fields of magnitudes 850 Oe and 450 
Oe. The results for H ext = 850 Oe arc multiplied by a 
factor of 10. As can be seen, the magnitude of this con- 
tribution is of the order of microvolts. It can exhibit a 
peak as a function of frequency. This peak is due to the 
fmr resonance of the permalloy metallic film when a' is 
close to zero. The magnitude of this peak depends on 
the Gilbert damping parameter which we have set to a 
value of 0.03. 

It is likely that the spins of YIG and that of the ferro- 
magnet are exchange coupled. The continuity of H x and 
By across the YIG-ferromagnet interface imposes partial 
alignment of the spins between these two materials. Spin 
currents in the FM will generate some residual spin wave 
in YIG as well. In the present calculation the exchange 
coupling between spins can be incorporated in the mag- 
netic susceptibility in eq. (1) by including the exchange 
energy in Jo- 
in conclusion we study the nonlocal spin and charge 
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current generation in a finite metallic element on the 
surface of magnetic insulators with the magnetic surface 
plasmon (MSP), which is absorbed and cannot be re- 
flected. The MSP generates a voltage in the longitudinal 
direction proportional to the power. If the metal is a 
ferromagnet, a spin current can also be generated in the 



longitudinal direction. Our result will be useful for new 
classes of devices using spin waves. 
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